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1 ‘^“SUMMARY. wDSO 


This is a computer program utilization report for the Automatic Data Manipula- 
tion Program (ADMAP), which is a single-precison FORTRAN IV aerospace data 
manipulation computer program developed by Analytical Mechanics Associates, Inc. , 
(AMA, INC.), to aid in the processing, reduction, plotting, and publication of electric 
propulsion trajectory data generated by the low thrust optimization program HII.TOP. 
ADMAP is set up to run on the IBM 360/91 at NASA Goddard. It has the option of 
generating SC4020 electronic plots, and therefore requires the SC4020 routines to be 
available at execution time (even if not used). It contains several general routines, 
including a cubic spline interpolation routine, electronic plotter dashed-line drawing 
routine, and single-parameter and double-parameter sorting routines. It also has 
many routines which are tailored for the manipulation and plotting of electric propulsion 
data, including an automatic scale -selection routine, an automatic curve-labelling 
routine, and an automatic graph-titling routine. It accepts data from either punched 
cards or magnetic tape. ADMAP is currently available in source-deck form to anyone, 
who might profit from its use. This report contains the basic instructions which 
enable a user to operate the program. 
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INTRODUCTION 

The ADMAP software package was created in 1969 to process the data which 
are output from the low-thrust trajectory optimization program HILTOP. The HILTOP 
program, which is described in Reference 4, came into existence during 1968-1969 as 
an improved version of a program developed originally at Princeton University. 
Because of HILTOP’s powerful ability to generate data, a relatively large amount of 
low-thrust trajectory and performance data began to come into existence during 1969 
at the Maryland office of Analytical Mechanics Associates, Inc. A need therefore 
arose to transform the data into a presentable and useful form, and ADMAP was con- 
ceived to fulfill this purpose. The final presentation format which was chosen for 
displaying low-thrust trajectory and performance data is exemplified by References 
1, 2, and 3. 

ADMAP is primarily designed for limited use by the author, to aid in the pub- 
lication of solar-electric propulsion data. As such, ADMAP experiences relatively 
little usage, and therefore a detailed documentation of the program is not warranted. 
The effort involved in preparing this user’s manual was therefore limited, and this 
report is essentially an unpolished first draft. 
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PROGRAM CAPABILITIES 


AD MAP's capabilities fall into two categories, those involving the SC4020 
electronic plotting routines and those not. 

When SC4020 plots are not desired, the program may be used to perform several 
basic data manipulation functions. The data, which appear in units called "data sets, " 
to be defined later, may be sorted with respect to any specified mission parameter, 
and redundant data will be eliminated or discarded. The data may be printed on the 
high-speed printer in two formats, one format being a straightforward listing of the 
data in condensed form and the other format being a display of the maxtmum and 
minimum values of each datum within a given data set. These features are available 
for data input to ADMAP by either cards or magnetic tape. The program may also 
generate a "tape catalog, " which consists of a brief summary of a magnetic data 
tape's contents which is output on the high-speed printer. A magnetic data tape is 
generated by ADMAP using card-data which have been punched by the trajectory op- 
timization program HILTOP. 

When SC4020 plots are desired, two plotting formats are available. The graphical 
format may be either in the general, one -graph-per -frame style or in the more specific, 
two -graphs -per -frame style as found in References 1, 2, and 3. The latter format 
is less general than the former since it was designed for the specific application of 
generating finished plots which are suitable for publication. The one -graph-per-frame 
format allows plotting of any specified subset of the data in either linear, logarithmic, 
or semi -logarithmic form, with automatic curve-labelling and graph-titling. In both 
formats, purely-balllstic missions are represented by dashed lines whereas low thrust 
missions are represented by solid lines, and all curves are generated by fitting a cubic 
spline function to the data. 


PROGRAM INPUT 


Inputs to ADMAP are of three basic types: 

(1) Card data-set input 

(2) NAMELIST input 

(3) Magnetic tape input 

Card Data Set Input 

A card "data set" is the basic unit of input to ADMAP. 

A card data set consists of 

(1) A title card 

(2) A NAMELIST input data group (defined later) 

(3) HILTOP trajectory summary cards 

(4) A special end-card consisting of five dollar-signs in the first five card 
columns and blanks elsewhere. 

The deck set-up of a card data set is as follows: 
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A card data-set thus contains NAMELIST input. A card data set is identical to a 
tape data set except that a tape data set contains no NAMELIST component. The HILTOP 
trajectory summary cards must be present in units of three cards or five cards per 
trajectory summary. At least one such trajectory-summary-unit must be present, 
and a maximum of 100 trajectory-summary-units are allowed per data set. The 
trajectory summaries throughout a given data set must all he comprised of either three 
or five cards each (No mixing is allowed.). The title card may contain any information 
whatsoever for the purpose of describing the given data set. The contents of the II It. TOP 
trajectory summary cards are given in Reference 4. The NAMELIST input data group, 
which must be present in each card data set, contains information which dictates to 
ADMAP how the data set is to be manipulated by the computer. 


NAMELIST Input 


Inputs to ADMAP are given through the NAMELIST feature of the IBM FORTRAN 
IV programming language. The input NAMELIST is named IN, and every control input 
required or used in the program is declared by name in the list. The general form 
for assigning an input value to a quantity is, simply 

NAME = VALUE 

where NAME is the FORTRAN name assigned to the quantity, which is included in 
the program's internal NAMELIST statement, and VALUE is a number consistent 
in form with NAME (i. e. , real, integer, or logical). All NAMELIST input cards 
must be blank in card column one. Unless otherwise specified, all NAMELIST names 
commencing with the letters I through N represent integers, whereas all names 
commencing with the letters A through H and O through Z represent floating-point 
numbers. All input data groups must begin with the characters &IN beginning in 
card column 2 and followed by a blank in card column 5, and end with the characters 
&END preceded by at least one blank. Multiple data assignments on a single card 
are separated by commas. The order of the input data assignments is arbitrary, and 
a quantity may be input more than once, the last value overriding all previous values. 
If no value is included in the inputs for a particular parameter, ADMAP s default 
value is used. For other details regarding the NAMELIST feature, see the IBM 
System/360 FORTRAN IV Language Manual. 

The following are the NAMELIST inputs. Values enclosed in parentheses denote 
default values. 

ITAPE Program master control flag. 

= (0) No magnetic tape is used. Numerical data is input by 

means of cards only. 

= 1 Hotting/listing of data which is input to the program 

by tape only. 

= 2 Generating a new data tape from data which is input 

by cards. 

= 3 Merging data from cards and old data tape to form a 

new data tape. 
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= 4 Copying an old data tape onto a new data tape and 

printing tape catalog. 

In the above description, "plotting" refers to electronic SC4020 plotting, and "listing" 
refers to the program's regular bond -paper printed output of data. This terminology 
will be used freely below. 

NTAPE Primary data tape logical unit assignment (17). 

NTAPE1 Secondary data tape logical unit assignment (18). 

When only one data tape is used, it is considered as the primary data tape. This 
tape is either being created as a new data tape from data cards, or is being used as 
data input to the program. When two data tapes are used simultaneously, the primary 
data tape is the "old" data tape which is only read, and the secondary data tape is a 
"new" data tape which is being written or generated. For example, NTAPE appears 
in the FORTRAN statement READ (NTAPE, 66) (OUT (I), 1= 1, 100), where 66 refers 
to a FORMAT statement and the OUT array consists of one trajectory -summary. 

LIST Plotting bypass flag. 

= (0) Hotting and listing of data. 

= 1 Listing of data (no plotting). 

LISTAL Generates listing of entire tape. 

= (0) Normal operating mode. Program is controlled by 

other inputs, notably NCODE, which is described below. 

= 1 Generates listing of entire (primary) data tape. Ignores 

plotting and other program inputs. 

ISORT Data sorting flag. 

= 0 Data is not sorted. H 

= (1) Data is sorted within each data set according to other .vj 

program inputs, which are described immediately j 

below. The amount of data constituting a "data set" is '(>, 

defined elsewhere. jl 

NSORT Primary-sort-parameter selection flag. Selects the parameter in ! 

the OUT array on which the current data set is to be sorted. Must 
have value 1, 2,3, ..... 100. Default value is 2. May be the 
negative of the above values only if the flag IPOW = 1, where IPOW 
is described below. The contents of the OUT array are given 
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TOL 


NSORT2 


TOL2 


explicitly in the description of the contents of an ADMAP data 
tape (sometimes referred to as a HILTOP data tape). The data 
in the data set, which consists of a collection of OUT-vectors, 
will be sorted according to ascending values of the NSORT-th 
parameter in the OUT-vector, that is, OUT (NSORT). The OUT- 
vectors which together comprise a data set are the objects which 
are actually sorted. Extra OUT-vectors having the same sort 
parameter value are deleted. (See TOL. ) 

Increment in value of primary sort parameter which determines 
when two sort parameter values are "equal" within this specified 
tolerance. Default value is . 099. After a data set has been sorted, 
extra OUT-vectors having the same value for the primary sort param- 
eter are deleted from the data set. This feature relieves the program 
user from worrying about possible duplicate data points in a given data 
set, since the program will automatically delete them. 

NOTE: A single OUT-vector contains a single HILTOP trajectory summary. 

Secondary-sort-parameter selection flag. Selects the parameter in the 
OUT array on which the current data set is to be sub-sorted (that is, 
sorted for the same value of the primary-sort-parameter). Default 
value is -10. Currently, 1POW must be equal to 1 in order for NSORT2 
to be used. NSORT2 is otherwise defined similarly to NSORT, discussed 
above. NSORT2 was designed for the purpose of sorting off-optimum- 
power data. Negative values of NSORT and NSORT2 cause sorting to 
be performed in descending rather than ascending order. 

Increment in value of secondary sort parameter which determines when 
two (secondary) sort parameter values are "equal" within this specified 
tolerance. Default value is .0099. TOL2 is defined similarly to TOL, 
discussed above, except that both primary and secondary sort parameters 
must be equal In two OUT-vectors (trajectory summaries) in order for 
one of the OUT-vectors to be eliminated from the data set. 


M 
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MIN MAX 


ITHREE 


LSCALE 


ICAT , 


NAMLST 


= (0) Normal operating mode. Program generates standard 
listing of data. 

= 1 Instead of standard listing of data, program generates a maximum- 
minimum map of the data, in which the maximum and minimum 
values of each parameter in a data set are displayed. 

(0) Each input trajectory-summary consists of five cards generated 
by the HILTOP program. 

= 1 Each Input trajectory-summary consists of three cards generated 
by the HILTOP program. 

ITHREE is used in subroutine READER only. 

Lagrange multiplier scaling flag. 

= (0) Do not scale the Lagrange multipliers. 

= 1 Scale the Lagrange multipliers such that the mass-ratio multiplier 
is unity. 

The scaling occurs immediately after reading either cards or tape. 

Tape catalog printing flag. 

= (0) Do not print tape catalog. 

= 1 Print tape catalog. A tape must be specified according to 

the input NTAPE. The Betting ITAPE = 4 is made internally. 

A tape catalog is merely a listing of title-records along with 
their code-numbers. 

NAMELIST suicide flag. - 

= (0) Read NAMELIST only once. 

= 1 Read NAMELIST until NAMLST = 0 occurs. 

NAMLST appears internally In the FORTRAN statement. 

IF (NAMLST, EQ. 1) READ (5, IN) 

NAMLST is to be used only when ITAPE=1 and the user desires to 

perform different functions with different tape data sets. It is easy 
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NCODE 


L1M1 


to get into trouble using NAMLST=1 and avoiding this is recommended. 

If NAMLST =1, each time a requested tape data set is encountered 
according to the input NCODE vector (described below), a NAMELIST 
input data group is also read into the program. The sequential NAMELIST 
data groups must correspond to the order of requested data sets as they 
appear on the tape, which is not necessarily the same order as specified 
in the NCODE vector. (A tape is read only once per job-submittal of 
ADMAP, with no backspacing over data sets.) 

Data-set code vector, dimensioned 500. Default values are all zero. 
NCODE is used only when ITAPE = 1, 2, or 3. When ITAPE - 2 or 3, 
input NCODE (1) = i in each NAMELIST input data group, where one 
NAMELIST input data group is contained in each card data set, and i 

z 

is a different (positive integer) code number for each data set. NCODE(l) 
therefore specifies the code number of the tape data set which is to be 
generated from the corresponding card data set. When ITAPE ~ 1. only 
those tape data sets specified by the NCODE vector are considered, all 
other tape data sets being ignored. For example, NCODE = 70, 12, 

1, 0, 4, 7314 would cause the manipulation (according to other NAMELIST 
inputs) of tape data sets having code numbers 70, 12, 1, 4, and 7314. 
Zeros are ignored; in fact, all zeros are eliminated by collapsing the 
NCODE vector such that no zeros appear between non-zero entries. 

The entries in the NCODE vector need not be in the same order as the 
sequential tape data sets, and NCODE-values not on the tape wifi cause 
no harm. A zero code-number indicates the end of the tape (which always 
consists of one file). When a tape is generated, a two-record data set 
having a code number zero is automatically generated at the end of the 
tape file. 

NCODE bypass (used only when 1TAPE=1). Default value is zero. When 
nonzero, L1M1 must be a positive integer used in conjunction with LIM2 
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LI M2 


IDOC 


IPUB 


IPOW 


JPOW 


(see below). When LIM1 is a positive integer, the NCODE vector is 
constructed internally to have values LIMl through LIM2 inclusive. 

This feature alleviates the necessity of inputing a large number of 
values into the NCODE vector. 

Upper code vector limit, as described immediately above. If LIM2 is 
less than LIMl, LIM2 is set equal to LIMl internally. Default value 
is zero. 

Electronic plotting master control flag. 

= (0) One-graph-per-frame plotting format. 

= 1 Two-graphs -per -frame plotting format, similar to the formats 

of References 1, 2, and 3. 

Electronic plotting master control flag (recently created). 

= (0) Normal program operation. 

= 3 Invoke special Initialization and logic for generating the electronic 
plots of Reference 3. EDOC is set equal to 1 internally. 

Control flag for manipulating special data sets which contain data for 
sub-optimal powered spacecraft. 

= (0) Data seta (card or tape) to be manipulated do not contain power- 
variation data. 

= 1 Data sets to be manipulated contain data in which the spacecraft 
reference power is the main variable. If IDOC=l, special logic 
for generating the electronic plots of Reference 2 is invoked. The 
two-parameter sort routtne currently cannot be used unless IPOW=l. 

Electronic plotter power normalization Indicator, used only when IDOC^l 
and IPOW =1. 

= (0) Plot % of optimum reference power on x-axis. This option assumes 
that the power has been sorted according to descending values of power 
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(for a given flight time), and that the first-power-value in the 
descending sequence is the optimum power. Reference 2 pro- 
vides an example. 

= 1 Plot actual reference power on x-axis, such as Reference 3 
illustrates. 

The following NAMELIST inputs apply exclusively to the control of the electronic 

plotting. 

IFIG Starting figure number when IDOC=l and sequential positive-integer 

figure numbers are desired. Default value is zero. First figure 
number will have the value IFIG +1. 

SCALE(I), I = 1,200 Curve scaling factors. Up to 200 curves may be plotted 

per frame when IDOC=0 and up to 200 curves may be plotted per two- 
frames when IDOC = 1. SCALE (1) is the scaling factor of the Ith 
plotted curve; that is, if SCALE (I) = 10, for example, then the Ith 
plotted curve has the real value obtained by multiplying the apparent, 
plotted value by 10. Default values are all 1. SCALE is ignored 
when LAUTO = 1. 

IAUTO Automatic scale factor computation flag. 

= (0) Curve scaling is performed by input curve scaling factors 
SCALE (I). 

= 1 Curve scaling Is performed internally by ADMAP in such a 

manner that each curve is plotted so as to fill out the frame 
(to the nearest whole order-of-magnitude). IAUTO is set to ( 

1 internally whenever IDOC = 1 . 

PC ENT Percentage of allowable overplotting (beyond the top and bottom graph 

borders) when IAUTO=l. Default value is 10. Portions of the curves 
lying above or below the borders are not plotted. 

LDASH Length (in raster units) of each dash for curve-segments in which 

the spacecraft reference power is zero. Default value is 8. 
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LGAP Length (in raster units) of each gap between dashes for curve-6egments 

in which the spacecraft reference power is zero. Default value is 4. 

NOBAL Ballistic-trajectory suppression flag for tape reading. ; 

= (0) Read all trajectory-summaries from an input tape. 

= 1 Read only those trajectory-summaries from tape into core- 

storage which are not ballistic (power >0). 

ICODEP Code number printing-flag. 

= (0) Do not print code number on each frame. 

= 1 Print data-set code number in bottom-left corner of each frame. 

1 t 

This option applies only when lDOC=l. 

The following NAMELIST inputs apply only when IDOC=0. ThiB option (1DOC -0) is 
not used very often, and some of the program logic under this option could be obsolete 
due to the several programming changes which have been made to ADMAP during the ]’ 

past year (1970) explicitly for the IDOC=l option. 

NPLQT The number of curves per frame. Default value is 1. 

NXAXIS X-axis parameter selection. Must be a positive integer. Default 

value is 2, which corresponds to flight time. The value of NXAXIS- i 
selects the l 4 * 1 element of the OUT-array as the quantity to be used as 
the independent variable (t. e. , x-axis variable) for plotting. The OUT- 
array is described In the section "ADMAP Data Tape Format and 
Contents. " 

* 

NYAXIS(I), 1= 1,200 Y-axis parameter selection. 

Default value of NYAX1S(1) is 8, corresponding to net spacecraft mass. 

NYAXIS(2) through NYAXIS (200) have no default values. The first 
NPLOT -elements of NYAXIS must be specified. NYAXIS is otherwise 
defined similarly to NXAXIS. 
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XMIN Minimum value on x-axis scale. Default value is zero. This quantity 

will appear at the leftmost point of the x-axis. 

XMAX Maximum value on x-axis scale. Default value is 8000. This quantity 

will appear at the rightmost point of the x-axis. 


DX Incremental value for x-axis scale. Vertical grid lines (if not suppressed) 

will appear along with their corresponding values along the x-axis at 
positions XMIN + DX, XMIN + 2DX, etc. Default value is 500. 


YMIN 

YMAX 

DY 


IT 

IB 

IL 

IR 


! 


NOLINE 


Defined similarly to XMIN, XMAX, and DX, respectively, except 
these quantities correspond to the y-axis. Default values are same as 
above . 

Raster locations of the top, bottom, left, and right graph borders, 
respectively. Default values are 960, 80, 100, and 980, respectively. 


Grid line suppression flag. 

= (0) plot grid lines. 

= 1 suppress grid lines and plot tic-marks instead. 


LINES The number of lines-deep which the letter-code section above each graph 

will be. Default value is 5. 


IBELOW Title-information position-adjustment flag. 

- (0) Print title Information just above the upper graph border. 

= 1 Print title information just below the upper graph border. 

The title -information comes from the title-card or title-record of the (first) data set 
from which information is being plotted. 


NOTELL Title-printing suppression flag. 
= (0) Print graph title. 

= 1 Do not print graph title. 
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NOLAB L Letter-code suppression flag. 

= (0) Print letter-code above each graph and next to each curve. 

= 1 Do not print letter-code. 

NWORD Graph title selection flag. 

= 0 No graph title is printed. 

= (1) "Solar -electric mission analysis" is printed. 

= 2 "Nuclear-electric mission analysis” is printed. 

NXCYL Logarithmic plotting flag for x-axis. 

= (0) Regular (non-logarithm ic) plotting 
= i , where i is a positive integer: 

x-axls will be logarithmic instead of linear, 
i is the number of base-10 log-cycles, and the smallest 
scale-value is determined by XMIN. Only whole-cycles are 
admissible. 

NYCYL Logarithmic plotting flag for y-axls, otherwise defined similarly 

to NXCYL. 

Note that plots may be log-log or semi-log. 

NSC The number of grid lines per log-cycle. Default value is 5. NSC 

applies to both x and y axes when logarithmic. 

SCALEQ(I), I = 1,10 Scale-values when plotting logarithmicly. The first NSC- 
elements of SCALEQ must specify the values (between 1 . and 10. ) of 
the cycle-subdivisions. Default values are 2. , 4. , 6. , 8. , and 10. , 
the remaining five elements being undefined. SCALEQ(NSC) must 
equal 10. for the graph to make sense. | 

MANY Data-set overplotting flag. Default value is 1. When MANY Is greater | 

than 1, NPLOT is set equal to 1, and the quantity specified by NYAX1S(1) ] 

will be overplotted from the next MANY data sets to produce a family of 
comparison curves of that quantity. £ 
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The above describes the inputs which may be used in each NAMELIST input 
data group, where by NAMELIST "input data group" is meant any card or set of 
cards commencing with &IN and terminating with &END. When using only a tape or 
tapes in conjunction with ADMAP (IT APE = 1 or 4), only one NAMELIST input data 
group is normally used, and this follows a single blank or dummy card. When card 
data sets are Input to ADMAP, a single NAMELIST input data group must be present 
immediately following the title card of each data set. (See the illustration describing 
the deck set-up of a card data set. ) 

The program is designed to have a minimum of NAMELIST inputs to perform 
a given task, since either the default values take care of the values not input or many 
NAMELIST inputs are computed internally as functions of other NAMELIST inputs. 
TheNAMELIST inputs used most often are ISORT, ITAPE, LIST, ICAT, IDOC, IPOW, 
LISTAL, and IPUB, and other inputs are seldom Invoked. 

Magnetic Tape Input 

Data may be input to ADMAP from a magnetic tape, and this option is indicated 
by using the NAMELIST input ITAPE. The contents and format of the tape must be as 
specified in the Program Output section. Usually, the input data tape has been previously 
generated from data cards using the ADMAP program with ITAPE=2. The data sets con- 
tained on the tape are manipulated according to the inputs specified in a single NAMELIST 
input data group which must be preceded by a single dummy-card. The dummy card is 
required as a substitute for a title-card of a (usually non-existent) card data set, since 
ADMAP thinks it is reading the first card -data-set until it encounters the first NAMELIST 
input data group and discovers that the input is by tape rather than cards. A £ape is read 
only once by the program, there being no capability at present to re-wind and re-read the 
tape according to subsequent NAMELIST specification. The program works with one data 
set (cards or tape) at a time, there being no capability at present to work with two or 
more data sets simultaneously. 
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Example Cases of Program Input 


The following page contains two examples of complete job input sequences for 
running ADMAP on the IBM 360/91 at the NASA Goddard Space Flight Center. 

The first input sequence uses no magnetic tapes but merely produces a standard 
condensed listing of the two input card data sets. Each data set is sorted acrording to 
ascending flight time by means of the default values of the NAMELIST inputs ISORT 
and NSORT. The input LIST=1 suppresses the electronic plotting and the input !THREE=1 
indicates that the HILTOP trajectory-summary cards are input as three cards per 
summary. The HILTOP trajectory-summary cards appear strange because they are 
coded in the A-format. 

The second input sequence provides an example of using ADMAP to generate 
electronic plots from data input by magnetic tape. The electronic plots are generated 
in the format of References 1, 2, and 3, as specified by IDOC=l, and only selected tape 
data sets are considered, as specified by the NCODE vector. Condensed listings are 
generated along with the plots. 


EXAMPLE CASE £1 


//yoifi JCE C^KC 

// EXEC LINKGC.REGICN.GC=340K 

//LI NK, SYSL IE DO O£N=£YS2.SC4O20 .C ISF=SFR 

/ /t tM< • S Y St i-fr-e e— * 

INCLUDE LCACL IE IN7F IMFLT ) 


ENTRY WAIN 

//GC .FT 1 IF CC 1 DC SX£CLT=F » CC E - (FECFM=FE .LRECL=8Q»ELK£ I2E=C0 ) 


//GO.OATA5 CC * 

1976 MCDE A CEPES FENC T A T ( 3 C ) / C EL 7 A / T E £ 6 A ( I A A C 1 ALFFA = £C 


-fc I N' t t S =E~1 t -M+yREir-i- -€ENe 

Y . 6 

) 

= • F 5 a) exl" 

< < - < 

4 

— -i ■ I- — -€>-■ rF- r > E — 


R 


Si 


/ 1 .A 


* 

|F .F 


C A , 


• i 

4 


C***3C7473C£CI 
0 £0 £ 12 32 

C 1C C73C£C£ 

S Q * * * 3 0 E 4 £ € C 2 C 1 

G 2C £ 12 32 

- ~€1 — 10 0^00303 

2 A G***209S63C£C1 

G 20 £ 12 32 

C 10 0 € £ C £0 3 


MODE A ENCKE PERI HE L I C N FENC T A T < 3 C ) / C EL T A / T E 2 6 4 ( 1 4 4 0 ) AtPFA = 30 


6 IK ti ST = i'« I THREE ISC'FHF*^ -£€-*€— 

A R A E 

, > 9 LC 


1 9C 


ieC*£€EC7CC £01 
160*3 6 €07 OX 20 2 
C60*3E6C70C££1 
O 40 £ 12 02 

C 1C C70C3C3 


- ***** 


EXAMPLE CASE #2 


-/-/YDL-R-0&6-GARC 

// EXEC L INKGC .REG l CN.GC = 340 K 

//LINK. SVSLIE DO C SN= S Y S2 • SG4O20 . C I SF-SFR 

//LINK.SY£L IN CD * 

INCLUDE LCACLI6 IN7E INFLT 1 
ENTRY MAIN 

- /-/■GO-^E-T-IGE 00 1---EE& — G-£-N-=-K7-R-|-Mt-I-T-A>N-I-T-— 2 4 0 3-7 .- - X 

// LABEL = ( .etF ) ,C ISF= (NEW .KEEP ) , V CLU W E = S E F= T C 77 7 C * X 

// DC6=( ELKS I 2E= 1 200 ,TRTCF = C.CEN= 1 > 

//GO .F T I 1FCCI CD £YSCLT = R >CCE- (F€CFM=FE . LRECL = EG» EL K£ I 2 E = EC ) 

//GO.FT l 7FOC 1 DO UN I T = 2A0C-9 . VGLL M£ = S £ fi=T C5 1 96; . L AE EL= < l.ELF ), X- ... 

// DI £P = ( CLC ,F ASS > ICCE = IRECFM = FE .ELKS I ZC=2 40 0 .CEN = 2 I 

'-•-/ /GG*DAT a»5— G C— * ; — 

DUMMY CARC 

61 N I OCC — 1 » NGC-CE = 7 , 3 8 « 2 *4 1 4 * 4 964 G 6 E N C 
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PROGRAM OUTPUT 


ADMAP outputs are of three basic types: 

(1) High-speed printer (unit 6) 

(2) Magnetic tape output (units variable). 

(3) Scratch pack output (unit 11) 

The high-speed printer output consists of either the standard condensed listing, a 
maximum-minimum map of parameters, or a tape catalog, as specified by various 
NAMELIST inputs. A tape catalog consists of a simple listing of the title-records of 
each tape-data-set together with its code-number. The magnetic tape output consists 
of a primary or secondary HILTOP data tape or an SC4020 plot -tape, as discussed in 
the NAMELIST inputs section. The SC4020 tape is on output unit 10. The primar\ 
data tape is on unit NTAPE and the secondary tape is on unit NTAPE1, where NTAPE 
and NTAPE1 are NAMELIST inputs. The scratch-pack output contains a brief job- 
summary. The contents of the high-speed printer condensed listing and the contents 
of an ADMAP (or HILTOP) data tape are given on the following pages. Examples of 
electronic-plot output may be found in References 1, 2, and 3. 
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Condensed Listing of Solar-Electric Mission Analysis Data 


Description of Format 


Header-descriptions which are not self-evident: 


PAY -LOAD 

V E H 
INIT MASS 
POWR PLNT 
PROP MASS 
TANK-AGE 
STRC TURE 
ALF 

POWER 1 AU 
PI, P2, P3 
PD1 , PD2, PD3 
PMASS 
ACCEL 
VJET 

V 001 

V 002 

CODE VECTOR 
RET PROP 
RET STRC 
RET ENGN 
ETA 

DEL VEL 
VEL LOSS 
PERI VEL 
HIGH BURN 


Net spacecraft mass at target planet, kilograms. 

Launch vehicle code-letter (may be ignored). 

Initial mass (launch vehicle payload), kilograms. 

Solar-electric propulsion system mass, kilograms. 

Solar-electric propellant ma6S, kilograms. 

Solar-electric propellant tankage mass, kilograms. 

Solar-electric structure mass, kilograms. 

Solar-electric propulsion system specific mass, kilograms/kilowatt. 
Reference power (at 1 A. U. from sun), kilowatts. 

Initial primer vector, Cartesian components, kg-tau/AU* 

Time derivative of Pi, P2, P3, kg/ A UJ * 

Mass-ratio initial Lagrange multiplier. 

Reference acceleration (at 1 A. U. from sun), meters/sec/sec. 
Solar-electric jet exhaust speed, meters/sec. 

Earth departure hyperbolic excess speed, meters/sec. 

Target planet arrival hyperbolic excess speed, meters/sec. 

(may be ignored). 

Retro-stage propellant mass, kilograms. 

Retro-stage engine + tankage + structure mass, kilograms. 
Retro-stage engine mass, kilograms. 

(may be Ignored). 

Retro-maneuver incremental velocity, meters/second 
Retro-maneuver velocity loss, meters/second 
Speed at perlapse of target planet capture orbit, meters/second 
Burn time of retro-stage, seconds. 


♦Coordinate system defined in Reference 3. 
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RET THR 

Retro-stage thrust, pounds. 

RET ISP 

Retro-stage specific impulse, seconds. 

R PER 

Capture orbit periapse distance, planet radii. 

R APO 

Capture orbit apoapse distance, planet radii. 

POWER TARG 

Solar-electric power at the target planet, kilowatts 

POWER MAX 

Maximum solar-electric power along trajectory, 
kilowatts. 

LOW BURN 

Low -thrust burn time, days. 

TRV ANG 

Total heliocentric travel angle, degrees. 

R MAX 

Trajectory maximum distance from sun, 
astronomical units. 

R MIN 

Trajectory minimum distance from sun, 
astronomical units. 

THR 

Solar-electric thrust at one astronomical unit 
from the sun, newtons. 

THR MAX 

Maximum solar-electric thrust along trajectory, 
newtons. 

R FINAL 

Spacecraft distance from sun at final time, 
astronomical units. 


Addendum - December 1970 

The above describes the data used in the publication of Reference 1. In actuality, 
the words "solar -electric" may be replaced by "low-thrust" in the above description 
of names, and the words "target planet"may be replaced by "final time" in many places. 
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ADMAP Data Tape Format and Contents 

An ADMAP data tape contains trajectory-summaries of converged optimum 
low thrust missions.* Each record consists of 100 words, and is created as follows: 
WRITE (NT A PE, 1) (OUT (I), 1 = 1, 100) 

1 FORMAT (100A4) 

Thus each word is of length 4 and is in the A-format, and the tape must be read the 
same way. A "data set" is defined to be a collection of records (which are always 
together and sorted according to ascending flight time) which describe a given mission 
(for example. Mode B Uranus orbiter missions with jettisoning of low thrust propul- 
sion system prior to capture), together with a title record, which is always the first 
record of the data set, the first 20 words of which contain the mission description 
(in shortened English). The contents of the title record are self-explanatory, except 
for orbiter missions, where "JET" means the low thrust system is jettisoned, and the 
absence of "JET" means It is not jettisoned. The last two records of the tape consist 
of words which are all zero. The OUT-array, mentioned in the WRITE statement 
above, is sufficient to describe each tape record (which contains data rather than title 
information). Letting (I) denote integer quantities and (R) denote real, or floating 
point, quantities, the OUT-array consists of the following:** 

OUT(l) = Julian launch date, days (xxxxx. x), (R). 

OUT (2) = Flight time, days (R). 

OUT (3) = Initial spacecraft mass, or payload of launch vehicle, kg (R). - 

OUT (4) = Propulsion system mass, kg (R). 

OUT (5) = Propellant mass, kg (R). 

OUT (6) = Propellant tankage mass, kg (R). 

OUT (7) = Structure mass, kg (R). 

OUT (8) = Net spacecraft mass, kg (R). 

OUT (9) = Propulsion system specific mass, kg/kw (R) 

OUT (10) = Reference power (at 1 AU from sun), kw (R). 

♦Generated by the HILTOP program, described in Reference 4. 

** Thus, for example, OUT(7) describes the 7th word of each tape record. 
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OUT(ll) = Power at target, kw (R). 


OCT (12) 
OUT (13) 
OUT (14) 
OUT (15) 
OUT (16) 
OUT (17) 
OUT (18) 
OUT (19) 
OUT (20) 
OUT (21) 
OUT (22) 
OUT (23) 
OUT (24) 
OUT (25) 
OUT (26) 
OUT (27) 
OUT (28) 
OUT (29) 
OUT (30) 
OUT (31) 
OUT (32) 
OUT (33) 
OUT (34) 
OUT (35) 
OUT (36) 
OUT (3 7) 
OUT (38) 
OUT (39) 
OUT (40) 
OUT (41) 


= Reference thrust (at 1 AU from sun), newtons (R). 

= Radial distance from sun at final time. AU (R). 

= Propulsion system efficiency (R). 

= Trajectory maximum distance from sun, AU (R). 

= Trajectory minimum distance from sun, AU (R). 

= MISHUN ( HILTOP flag which is now obsolete). (I). 

= MBOOST ( HILTOP launch vehicle flag), (I). 

= Maximum power along trajectory, kw (R). 

= Maximum thrust along trajectory, newtons (R). 

= Low thrust burn time, days (R). 

= Total heliocentric travel angle, degrees (R). 

= ETA (HILTOP efficiency flag), (R). 

= Retro stage engine mass, kg (R). 

= Retro stage engine + tankage + structure mass, kg (R). 
= Retro stage propellant mass, kg (R). 

= Burn time of retro stage, seconds (R). 

= Retro stage thrust, pounds (R). 

= Retro stage specific impulse, seconds (R). 

= Capture orbit perlapse distance, planet radii (R). 

= Capture orbit apoapse distance, planet radii (R). 

= Capture orbit periapse speed, meters/sec (R). 

= Retro maneuver incremental speed, meters 'sec (R). 

= Retro maneuver velocity loss, meters/sec (R). 


Initial primer vector, Cartesian components.* 
kg-tau/AU (R). 


Time derivative of initial primer vector,* 
Cartesian components, kg/AU (R). 


= Mass ratio initial Lagrange multiplier (R). 


♦Coordinate system defined In Reference 3. 
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OUT (42) 
OUT (43) 
OUT (44) 
OUT (45) 
OUT (46) 
OUT (4 7) 
OUT (48) 
OUT (49) 
OUT (50) 
OUT (51) 
OUT (52) 
OUT(53) 
OUT (54) 
OUT (55) 
OUT (56) 
OUT (57) 
OUT (58) 
OUT (59) 
OUT (60) 
OUT (61) 
OUT (62) 
OUT (63) 
OUT (64) 
OUT (65) 
OUT (66» 
OUT (67) 
OUT (68) 
OUT (69) 
OUT (70) 
OUT (71) 


- Reference acceleration (at 1 AU from sun), meters/sec/sec (R). 
= Jet exhaust speed, meters/sec (R). 

= Earth departure hyperbolic excess speed, meters/sec (R). 

= Initial time, days (R). 

= Final time, days (R). 

= Target planet arrival hyperbolic excess speed, meters/sec (R). 

= blank. 

= blank. 

= blank. 

= MODE ( HILTOP power variation flag), (I). 

= MOPT (HILTOP first guess option flag), (I). 

= MOPT2 (HILTOP initial planet flag), (I). 

- MOPT3 (HILTOP target planet flag), (I). 

= MTMASS (HILTOP retro propulsion flag), (I). 

= IEPHEM ( HILTOP analytic ephemeris flag), (I). 

= JPP (HILTOP propulsion system jettison flag), (1). 

= JT (HILTOP tankage jettison flag), (I). 

0 

= V Final primer vector, Cartesian components, kg-tau/AU (R). 


Time derivative of final primer vector, Cartesian 
components, kg/AU (R). 


= Mass ratio final Lagrange multiplier (R). 


Final spacecraft heliocentric position, Cartesian 
components, AU (R). 



Final spacecraft heliocentric velocity, Cartesian 
components, AU/tau (R). 


23 


OUT(72) = Final mass ratio (R). 
OUT (73) = 


Launch planet position at initial time, Cartesian 
components, AU (R). 


Launch planet velocity at initial time, Cartesian 
components, AU/tau (R). 


OUT(74) = 

OUT (75) = 

OUT (76) = 

OUT (7 7) = 

OUT (78) - 

OUT (79) = Initial time, tau (R). 

OUT (80) = 

OUT (81) = 

OUT (82) - 
OUT (83) = 

OUT (84) = 

OUT (85) = 

OUT (86) = Final time, tau (R). 

OUT (87) through OUT(99) = blank. 

OUT (100) = data set or mission code number (1). 

The integer code number is the last (100th) word of each record, including the 
title record, and is the same code number given in the TAPE CATALOG. A zero code 
number indicates the end of the tape. Each data set has a unique code number associated 
with it. 


Target planet position at final time, Cartesian 
components, AU (R). 


Target planet velocity at final time, Cartesian 
components, AU/tau (R). 


Addendum - December 1970 

The above describes the ADMAP data tape used in the publication of Reference 1. 
Since that time, two more reports (References 2 and 3) have been published using ADMAP, 
and their data tapes differ slightly from the Reference 1 data tape. Specifically, the 
latter tapes contain trajectory summaries describing sub-optimal powered spacecraft, 
and each data set (collection of tape records) which contains sub-optimal data has a title- 
record containing the mnemonic "POW". The latter tapes are blank (zero) in OUT(59) 
through OUT(86). All of the above comments concerning the title-record of a tape-data- 
set apply to the three data tapes in existence at this writing, those used to generate 
References 1, 2, and 3. In actuality, the title-record may contain any information 
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(80 characters) which the program user desires for the purpose of describing the 
information in the data set of which the title-record is a part. The tape -title-record 
is identical in content to the title-card of a card-data-set, as the title card is copied 
"directly" onto the tape at the time the tape is created from cards. 

The contents of the OUT-array, as given above, apply to all discussions through- 
out this report which make reference to the OUT-array. The OUT-array is the basic 
vector of quantities used for all purposes within the ADMAP program. 
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BRIEF SUBROUTINE DESCRIPTIONS 


Identification 


MAIN 


ADJUST 


AUTO 

CYCLE 


DASHV 


DASHW 


DISECT 


ENDCAS 


Purpose 

Program master control. Loads data from tape or 
cards into core array by calling READTP or READER. 
Sorts data by calling SORT or SORT2 and prints con- 
densed listing of data by calling PRINT. Supervises 
electronic plotting within a set of nested D‘ '-loops 
(statements 1000. 29, and 11). Terminates program 
execution. 

Adjusts the position of the x-axis scale numbers 
on the electronic plots so that they are centered 
with respect to the lines which they represent. 

Thus larger numbers do not appear lopsided. 

Automatically selects the scale-factors for the 
plotted functions. Called from MAIN program only. 

Reads tape data in which the user is not interested. 
Reads or cycles through all tape records having the 
same code number (100th word of record) until code 
number changes, then backspaces to start of next 
data set (which is defined by a code number)*. 

Provides plotting capability of drawing dashed lines. 
See descriptive comments at beginning of subroutine 
listing. 

Provides capability of specifying horizontal and 
vertical boundaries for curves plotted using subroutine 
DASHV. 

Decomposes a floating-point number by shifting the 
decimal point until there is only one nonzero digit 
to the left of the decimal point, and computes the 
corresponding exponent. 

When reading data from cards, this routine tests for 
the special card ($$$$$) which flags the end of a 
"data set" (which is defined by the comment card 
preceding it)*. 


♦The terminology "data set" is defined in the PROGRAM INPUT section. 
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EXTRAP 

FINISH 

FUNCT 

GRID 

IDENT 

INPUT 

LABEL 

LINEW 

MAXMIN 


Performs linear or quadratic extrapolation of 
plotted functions to solve for the point where the 
reference power goes to zero (if such is the case). 

Finishes merging tape and cards to generate a 
new data tape, by copying the old tape directly 
onto the new tape but Ignoring data sets corres- 
ponding to the cards, which have been read onto 
the first part of the new tape. 

Converts the function to be plotted into logarithmic 
form when log-log or semi-log plots are requested. 

Generates the coordinate grid for each graph. 
Computes scaling factors between floating point 
and SC 4020 raster values of the quantities to 
be plotted. Computes SC4020 raster location of 
the origin of coordinates for each graph. 

Generates SC4020 identification frame when 
IDOC=0. 

Supervises program input by NAMELIST. Initiates 
SC4020 plotting. Controls program capability of 
overplotting functions from different data sets 
when JDOC-O. Initiates magnetic tapes. Sets 
default values of input quantities which are func- 
tions of other input quantities. 

Automatically labels plotted curves with SC4020 
alphabetic characters according to a preset code. 
Contains limited memory to avoid overplotting of 
labels. 

Provides capability of specifying horizontal and 
vertical boundaries for curves plotted using SC4020 
subroutine LrNEV. 

This routine is similar to subroutine PRINT, except 
it generates a map of the maxima and minima of 
the quantities in a given data set rather than a con- 
densed listing of data. 
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MOMENT 

PLOT 

POW 

PRCAT 

PRINT 

READER 

READTP 

SEARCH 


Computes moments for the cubic spline function. 

Supervises plotting of an SC4020 graph. Selects 
number of straight line segments which compose 
each curve. Prevents a flat curve from grazing a 
border. Determines segmentation of curve into 
solid and dashed components. 

Examines set of spacecraft reference power 
values to determine which subsets are zero, and 
supervises the extrapolation of function values to 
where reference power becomes zero (from the 
nonzero side). 

Prints a so-called "tape catalog" which consists 
of merely a listing of all of the title cards 
which precede each set of data, which is equivalent 
to the first tape record of each data set. This 
gives the user an idea of a data tape'B contents 
and the order in which the data appears on the 
tape. 

Generates the condensed listing of data (standard 
program printout). 

Reads five (or three) data cards each time called. 
The five (or three* cards represent a single HILTOP 
trajectory summary. Scales the Lagrange multi- 
pliers to unit mass multiplier if requested. 

Basic purpose is to transfer data from cards to 
core storage. 


Basic purpose is to read requested tape data into 
core storage, bypassing undestred data. Copies 
old data tape onto new one if requested. 

Searches the input set of code numbers NCODE to 
Bee If the current code number on the data tape 
(which is being read) matches any element of 
NCODE. If a match occurs, a flag (NFLAG) is 
set to 1. 
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SETUP 


SORT 


SORT 2 


SPLINE 


SPLIT 


START 


When IDOC~l, this routine sets up (computes) the 
various internal plotting parameters required for 
each of the four plots associated with a single 
figure. (IDOC=l triggers the logic which creates 
the plotting format used in the solar-electric pro- 
pulsion documents which have been published. ) 

Sorts an array so that Nth element of array is 
monotonically increasing. Eliminates all dupli- 
cate columns of the array by condensing the 
array. Sorting method is not efficient for very 
large arrays. 

Two-parameter sort routine. Sorts an array 
according to ascending or descending values of 
one element, then sorts each subset of the array 
(for which the first sort parameter has the same 
value) according to ascending or descending 
values of a second element of the array. Elimi- 
nates all duplicate columns of the array by con- 
densing the array. Sorting method is not efficient 
for very large arrays. 

Cubic spline function. This is the curve-fit used 
for each function which is plotted. For each 
function which is plotted, SPLINE must be initialized 
by a single call to subroutine MOMENT, which 
computes certain constants, called moments, which 
are associated with the spline curve. 


When a curve consists of more than one segment, 
which occurs when the spacecraft reference power 
vanishes, this routine sets up the x-axis and y-axis 
values which are to be plotted on the next curve - 
segment. Each curve-segment is plotted as though 
It were a lone curve, unassociated with other segments. 

This routine sets input quantities to their default 
values. Called once at start of execution only. 

This routine is equivalent to BLOCK DATA except 
it will run on any machine which accepts standard 
FORTRAN. 


29 



Addendum 


TAPE 


TITLE 


TIT LEX 


ZAP 


Generates sequential tape records. Basic purpose 
is to transfer data from core storage to magnetic 
tape. 

Prints (on SC4020 output tape) graph titles and 
curve label codes when the input flag IDOC O (One 
graph per SC4020 frame mode of operation). 

Prints (on SC4020 output tape) graph titles and 
curve label codes when the input flag IDOC’=l (Two 
graphs per SC4020 frame mode of operation). 

Decomposes a floating point number into a vector 
of hollerith characters so that the number can be 
printed using the SC4020 subroutine PRINTV. 


December 1970 

FIGURE Prints (on SC4020 output tape) the figure number 

of each figure. 

SETUP3 When IPUB=3, this routine sets up (computes) 

certain parameters required in the generation of 
the electronic plots for Reference 3. 
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MACHINE REQUIREMENTS 


AJDMAP's core requirements on the NASA Goddard IBM 360/91 computer are 
3E758 hexadecimal locations. The source program is single-precision USA standard 
FORTRAN IV and consists of about 3,500 cards. ADMAP requires that the standard 
library routines (ABS, ALOGIO, etc. ) be available and also that the SC4020 electronic 
plotter routines be available at execution time. If no electronic plotting is desired, 
the SC4020 routines could be replaced by dummy routines. 

For peripheral equipment, ADMAP uses the standard UNIT 6 for high-speed 
printed output, and UNIT 5 for input by the card reader. In addition to these, the 
program uses UNIT 11 for output onto the NASA Goddard Remote Input Terminal 
System (RITS). The UNIT 11 output consists of a brief job summary and may be 
diverted to a scratch pack if desired. The program uses a maximum of two magnetic 
tapes at any given time. At NASA Goddard, the SC4020 plot tape is required to be 
mounted on UNIT 10, although no unit is explicitly specified within ADMAP. If one or 
two HILTOP data tapes are to be manipulated, their tape units are specified by the 
NAMELIST entries NTAPE and NTAPE1. Specific examples of IBM job control cards 
for running ADMAP are found in the section "Example Cases of Program Input. " 

The majority of ADMAP runs are less than one minute of CPU time. The 
longest duration ADMAP runs to date have been the three runs which generated the 
electronic plots for References 1, 2, and 3, and these runs consumed from four to 
seven minutes of CPU time with the I/O time being somewhat less. 
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